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"<n": 

, We study rigidity transitions in covalent amorphous networks using size- 

. increasing cluster approximations and constraint counting algorithms. Pos- 

. sible consequences of the presence of self-organization are examined. The 

1-^ ' analysis reveals two transitions instead of the usual (mean-field) single tran- 

sition. One from a floppy to an isostatic rigid phase at a mean coordination 
number fd where the number of floppy modes vanishes and a second one from 
^ \ an isostatic to a stressed rigid phase at fc2 . The value of the two critical mean 

, coordination numbers as well as the width Af = fc2 — ^ci of the intermediate 

. phase depend very strongly on the presence of medium range order elements 

O ■ such as rings. 

Q ; Pacs: 61.43.Fs-46.30.Cn 

e: 

\ The notion of constraints and their application to classical macroscopic physics problems 
O ' such as the stability of bridges and trusses have been introduced and first considered by J.L. 
^ . Lagrange and J.C. Maxwell On this basis, J.C. Phillips asserted some twenty 

' years ago that covalent networks can be mechanically constrained by interatomic valence 
' forces such as bond- stretching and bond-bending and optimal glass formation is attained 
■ " " ' when the netwok sits at a mechanically critical point. This happens when the constraints 
Tic per atom estimated by Maxwell counting equal the degrees of freedom per atom in 3D, 
i.e. ric = 3. 

Such mechanical systems have been examined in terms of percolation theory by M.F. 
Thorpe [^] who showed by a normal mode analysis that the number of zero frequency solu- 
tions (floppy modes) / of the dynamical matrix equals f = 3 — Uc and vanishes when the 
mean coordination number f of the network reaches the critical value = 2.4. In this mean- 
field approach, one considers a network of N atoms composed of rir atoms that are r-fold 
coordinated. The enumeration of mechanical constraints in this system gives r/2 bond- 
stretching constraints and (2r — 3) bond-bending constraints for a r-fold coordinated atom. 
Since then, a certain number of structural possibilities have been taken into account such as 
rings, broken bond-bending constraints or the effects of one-fold coordinated atoms [0. 
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These powerful ideas have led to the prediction of a floppy to rigid transition in random net- 
works and various examples where rigidity percolation threshold occurs have been reported 
1^. Also, applications of rigidity in biology and computational science have been reported 
H], |10|. Nevertheless, experiments on binary and ternary chalcogenide glasses have shown 



the existence of two transitions at Vd and instead of the single mean- field transition |11 



P^ . This suggests that the mean-field constraint counting alone as it has been realized up 
to now, may be insufficient to describe accurately the underlying phase transitions. Recent 
results on Raman scattering and modulated differential scanning calorimetry (MDSC) real- 
ized on Group IV binary chalcogenides or ternary glasses suggest indeed evidence [0] for 
the growth of a self-organized (isostatic rigid) intermediate phase between the floppy and 
the stressed rigid phases, for which evidence is obtained from numerical simulations |jl5| . 

However, a certain number of questions remain at this stage. What controls the values fd 
and fc2, the width Af = fc2 — rd ? Recent results show that this width can be particularly 
sharp |jl6| . How does isostatic regions and self-organization influence the absolute magnitude 
of these quantities ? What can be done which goes beyond the elegant mean-field approach 
? This Letter attempts to adress these basic issues. We report here on the role of medium 
range order (MRO) in glasses of the form B^Ai^x with coordination numbers = 2 and 
rs = 4, and f = 2 + 2x. Typical glasses are the Group IV chalcogenides such as GexSei^^ 
which have been extensively studied in this context. To construct MRO, we have used 
size-increasing cluster approximations (SICA) to generate sets of clusters on which we have 
applied constraint counting algorithms. The results show two transitions, one at which the 
number of floppy modes vanishes. Another transition (a "stress transition") where stress 
in the structure can not be avoided anymore, is located beyond. In between, this provides 
evidence for a self-organized network for which the probability of stress-free clusters has 
been computed. The width Af increases with the fraction of MRO elements. Finally, in 
case of random bonding, a single transition is obtained. 

Construction. SICA have been first introduced to elucidate the formation of fuUerenes 



17|, but also the intermediate range order in amorphous semi-conductors [IS]. They rely 



on the statement that the fraction of significant MRO structures converges very rapidly to 



a limit value when the size of the considered clusters is increasing |]T9|. The construction is 
realized in Canonical Ensemble with particular energy levels. One starts from short range 
order molecules (the basic units at the initial step / = 1 which will serve as building blocks) 
and construct all possible structural arrangements of two basic units (/ = 2, see Table I), 
three basic units (/ = 3) and so on. This is supposed to be realized at the formation of 
the network, when T equals the fictive temperature Tf [^. Here, we have chosen as basic 
units the A2 and the stoichiometric balanced BA2 molecules (e.g. Se2 and GeSe2) for a 
reason which will become clear below. We have checked that the results do not depend 
on this particular initial choice. These basic units have respective probabilities 1 — p and 
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j9 = 2x/(l — x), X being the concentration of B atoms. 

The creation of a chain-hke A2 — A2 structure will involve an energy gain of Ei, isostatic 
A2 — BA2 bondings will use an energy gain of E2 and the creation of corner-sharing (CS) and 
edge-sharing (ES) -8^4/2 tetrahedra respectively £'3 and £'4. The latter quantity will be used 
to tune the fraction of ES among the structure. The produced probabilities have different 
statistical weights which corresponds to the number of equivalent ways a given cluster can be 
constructed. This quantity can be regarded as the degeneracy of the corresponding energy 
level. For instance, given the coordination number 4 of the basic unit -8^4/2 and labeled 
covalent bonds a CS -82^4 cluster has the multiplicity 4x4, whereas for a ES cluster, we 

count 2 X = 72 in three dimensions 

Due to the initial choice of the basic units, the value of the energy E2 will influence 
the probability of isostatic clusters since this quantity is involved in the probability of the 
isostatic BA4 cluster {ric = 3, see Table I). If we have 62 ^ 61,63,64, the network will be 
mainly isostatic. 

At step / = 2, we can generate three types of clusters (Table I), A4, BA^ and i?2^4 having 
two isomers (the CS and ES tetrahedra). Their unrenormalized probabilities are given by: 
PA^ = 4(1 - p)^6i, PBA4 = 16p(l - p)e2. Pes = 16p^63 and pss = 72^^64 out of which can 
be extracted the concentration x^'^^^ of B atoms. The quantities 6j = exp[—Ei/Tf] are the 
Gibbs weights at Tf. Next, we compute the number of mechanical constraints (bond-bending 
and bond-stretching) per atom on each cluster by Maxwell counting. Special care has to be 
taken in order to avoid the counting of redundant constraints on clusters containing rings, 
following the procedure described by Thorpe The probabilites depend on two parameters 
(i.e. the Gibbs weights ei/62 and 63/62) and eventually 64/62 if one considers the possibility 
of ES or rings. One of these two weights can be calculated by writing a conservation law for 



the concentration of B atoms [21 



X 



(') = x (1.11 



These weights become composition dependent in solving equ. ( [L.l|) which means that either 
the energies Ei or the fictive temperature Tf depend on x but here only the 6j(x) 
dependence is relevant for our purpose. With increasing cluster size, it is obvious that the 
number of potential isomers will increase (Table I), also the different types of rings which 
have some evidence in chalcogenides We have realized the construction up to the step 
1 = 4. At each step, we have determined either 61/62 or 63/62 solving equ. (|1.1| ) and 



computed the total number of constraints ric per atom on the set of clusters (see Table I). 
Finally, we have looked for the concentration of B atoms (or the mean coordination number 
f) for which the number of floppy modes vanishes. 

Results. Random bonding is obtained by setting the above defined Gibbs weights 6^ to 
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one and the cluster probabilities are then only given by their statistical weights. 

Solving / = 0, one obtains a single transition for all steps in the mean coordination 
number range [2.231,2.275], somewhat lower than the usual mean-field value 2.385 |TH]. This 
comes from the fact that the number of equivalent ways to connect units together 

is substantially higher than for the connection of (chain-) A2 units. We do not obtain an 
intermediate phase for random bonding. 

Let us turn to self-organization and proceed as follows. Starting from a fioppy cluster 
of size / (almost a chain-like structure made of A atoms), we allow the agglomeration of a 
new basic unit onto this cluster to generate the cluster of size / + 1 only if the creation of 
a stressed rigid region can be avoided on this new cluster (due to the agglomeration of a 
BA4/2 basic unit onto another BA4/2 tetrahedron being part of the /-sized cluster). With this 
rather simple rule, upon increasing f we will accumulate isostatic rigid regions on the size 
increasing clusters because -8^4/2 units are only accepted in A2 — -BA4/2 isostatic bondings. 
Alternatively, we can start from a rigid cluster which exist at higher mean coordination 
number (f < 2.67) and follow the same procedure but in opposite way, i.e. we allow only 
bondings which lead to isostatic rigid regions, excluding systematically the possibility of 
fioppy A2 — A2 bondings. 

Here, the simplest case deals with dendritic clusters, where we have removed all possi- 
bilities of ring creation. For / — > 00, this would permit to recover the results on Random 
Bond Models |^3|] for which there are no loops or rings in the thermodynamic limit and to 
obtain equivalence with Bethe lattice solutions [^|. We obtain a single transition for even 
/ steps at exactly the mean-field value f = 2.4 whereas for the step / = 3, there is a sharp 
intermediate phase defined by / = (again at f = 2.4) and the vanishing of fioppy regions 
(i.e. ei/e2 is zero) at f = 2.382(6). The probability of isostatic clusters as a function of the 
mean coordination number has been computed and shows that the network is entirely stress 
free at the point where / = (solid line, fig. 2). If there is a width (for / = 3), then the 
same probability is less than one and displays a narrow distribution. 

Next, we have allowed a certain amount of medium range order (MRO) by setting the 
quantity 64/62 7^ 0. Two transitions are then obtained for every SICA step. A first one at 
fci where the number of fioppy modes vanishes. A second one at fc2 defined by 63/62 = 0. 
This means that beyond this point, stressed rigid regions created by the connection of at 
least two BA4^/2 units can not be avoided anymore, i.e. the Gibbs weight 63/62 becomes 
non-zero, and composition dependent. We call this point the "stress transition" because its 
definition is very close to the one given in [15]. We show the / = 2 result (fig. 2) where 
/ = at fci = 2.4 and 63 7^ at f > fc2 for different fractions of ES tetrahedra, defining 
the intermediate phase Af. f^i does not depend on the ES fraction, as well as the fraction 
of stressed rigid clusters in the structure. To ensure continuous deformation of the network 
when B atoms are added and keeping the sum of the probability of fioppy, isostatic rigid 
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and stressed rigid clusters equal to one, the probability of isostatic rigid clusters connects 
the isostatic solid line at fc2- Stressed rigid rings first appear in the region fd < f < fc2 
while chain-like stressed clusters (whose probability is proportional to 63) occur only beyond 
the stress transition, when 63 7^ 0. We conclude that when f is increased, stressed rigidity 
nucleates through the network starting from rings. Results remain similar for the even 
/ = 4 step. It appears from fig. 2 that the width Af = fc2 — fci of the intermediate phase 
increases with the fraction of MRO. We have represented this quantity as a function of the 
MRO fraction at the rigidity transition in fig. 3 which shows that Af is almost an increasing 
function of the ES fraction as seen from the result at SICA step / = 4. Here, there is only 
a small difference between allowing only four-membered rings (ES) (lower dotted line) or 
rings of all sizes (upper dotted line) in the clusters. Finally, one can see from fig. 2 and the 
insert of fig. 3 that the probability of isostatic clusters is maximum in the window Af, and 
almost equal to 1 for the even SICA steps, providing evidence that the structure is almost 
stress-free. 

Discussion. Chalcogenide glasses represent the ideal systems to check these results. 
Different types of experimental measurements have given evidence on the two transitions 
and the nature of the self-organized intermediate phase. Raman scattering has been used 



ri|,p!3[] as a probe to detect elastic thresholds in Si^Sei^x and Ge^Sei^x glasses. Specifically, 



changes in the CS mode chain frequencies have been studied with glass compositions and 
show a kink (or a jump) at the mean coordination number fd = 2.4 and fc2 = 2.52 in Ge and 
fc2 = 2.54 in Si based systems, suggesting onset of a new rigidity at fc2- A clear correlation 
between these results and the vanishing of the non-reversing heat fiow AHnr (the part of 
the heat fiow which is thermal history sensitive) in MDSC measurements has been shown 
|]Tl|,|T3| . Obviously, since this AHnr term provides a measure of how different a glass is from 
a liquid in a configurational sense, this suggest that in the intermediate phase, glass and 
liquid structure are closely similar to each other. 

The SICA and constraint counting algorithms show that the width Af of the intermediate 
phase increase with the fraction of ES tetrahedra and more generally with MRO composed 
of small rings (fig. 3). We stress that the width should converge to a lower limit value of Af 
compared to the step 1 = 2, therefore one can observe the shift downwards when increasing 
/ from 2 to 4. This limit value is in principle attained for / — > 00, or at least for much larger 



steps than I = 4 [T^. For Si-Se, Af = 0.14 is somewhat larger than forCe-Se (Af = 0.12) 
consistently with the fact that the number of ES is higher in the former [|1^ . 

In summary, we have shown that size increasing cluster approximations could be used 
to go beyond the mean-field approach of the rigidity transitions. We have estimated for 
the different approximation steps the number of mechanical constraints and the number of 
floppy modes / and shown that two transitions were obtained in this situation. One at 
which / vanishes and another at which stressed rigid regions appear on the clusters. The 
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width Af is an increasing function of the MRO fraction. In the window Af, the rate of 
isostatic clusters is at its maximum. These new results should motivate developments on 
the role of local structure and MRO in the rigidity transition, and applications to Group V 
chalcogenides such as As^Sei^x glasses. 

The Laboratoire de Physique Theorique des Liquides is Unite Mixte de Recherche n. 
7600 du CNRS. 
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3 Width of the transition Af as a function of the fraction of MRO clusters at 
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(dotted lines). At step 1 = 2, the MRO clusters reduce to the edge-sharing 
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having only ES as MRO element. 'I'he insert shows the probability of isostatic 
clusters with mean coordination number f for / = 4 (dotted line) and / = 3 
(dashed line). 'I'he shaded region of / = 4 is deflned by the corresponding Af" 

List of Tables 

I Clusters generated at the different SICA steps / with the chemical formula in 

case of GcxSei-x glasses, the number of isomers and the number of constraints 

He per atom. 'I'he number of clusters containing rings is indicated in bracketts. 
GeSe^ and Ge2Se^ are isostatic clusters with respective energy levels E2 and 
2^2.1 



8 




FIG. 1. Some of the Ge^Se^ and Ge2SeQ MRO clusters generated by SIC A at step I = 3 with 
their statistical weight, a) A six-membered ring with ric = 3.67 b) An edge-sharing tetrahedra 
chain typical of vitreous SiSe2 with ric = 3.22 c) A six-membered ring with chalcogen inclusions 
and ric = 3.25. 
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Mean coordination number r 

FIG. 2. Probability of floppy, isostatic rigid and stressed rigid clusters as a function of the mean 

coordination number for different fractions of ES at / = 2. The solid lines correspond to the 
dendritic case where no edge-sharing tetrahedra are allowed. The broken lines correspond to the 
same quantities for ES fraction at the stress transition of 0.156, 0.290 and 0.818. For a ES fraction 
of 0.156, the filled squares indicate the point fc2 at which the stress transition occurs and serves 
to define the intermediate phase Af. 
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Fraction of MRO at the rigidity transition 

FIG. 3. Width of the transition Af as a function of the fraction of MRO clusters at the rigidity 
transition for Z = 2 (sohd hne), I = 3 (dashed Une) and I = 4 (dotted Unes). At step I = 2, the 
MRO clusters reduce to the edge-sharing GeS 64^/2 tetrahedra. For larger steps, different rings sizes 
(4,6,8) have been taken into account. The lower dotted line correspond to a system at Z = 4 
having only ES as MRO element. The insert shows the probability of isostatic clusters with mean 
coordination number f for Z = 4 (dotted line) and I = 3 (dashed line). The shaded region of Z = 4 
is defined by the corresponding Af. 
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Size I 


cluster 


Number of isomers 


ric 


1 


Se2 


1 


2 




GeSe2 


1 


3.67 


2 


Se4 


1 


2 




GeSe^ 


1 


3 




Gt^Sci 


2(1) 


3.67 


3 


See 


1 


2 




GeSee 


2 


2.71 




Ge2Se6 


4(2) 


3.25 




Ge^See 


4(3) 


3.67 


4 


Ses 


1 


2 




GeSes 


3 


2.56 




Ge2Ses 


11(6) 


3 




G 63 Ses 


12(9) 


3.36 




Ge^Ses 


10(9) 


3.67 



TABLE I. Clusters generated at the different SICA steps / with the chemical formula in case of 
GcxSex-x glasses, the number of isomers and the number of constraints ric per atom. The number 
of clusters containing rings is indicated in bracketts. GeSe^ and Ge2Ses are isostatic clusters with 
respective energy levels E2 and 2E2. 
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